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Studying various thermodynamic quantities for the free domain waU fermions for both finite and 
infinite fifth dimensional extent N5, we find that the lattice corrections are minimum for Nt > 10 
for both energy density and susceptibility, for its irrelevant parameter M in the range 1.45-1.50. The 
correction terms are, however, quite large for small lattice sizes of Nt < 8. We propose modifications 
of the domain wall operator, as well as the overlap operator, to reduce the finite cut-off effects to 
within 10% of the continuum results of the thermodynamic quantities for the currently used Nt ~6- 
8 lattices. Incorporating chemical potential, we show that divergences are absent for a large class 
of such domain wall fermion actions although the chiral symmetry is broken for fJ, ^ 0. 

PACS numbers: ll.15.Ha, 12.38.Mh, 12.38.Gc 



I. INTRODUCTION 

The nature of QCD matter at high temperatures and 
densities has been of interest due to the experiments at 
the Relativistic Heavy Ion Collider(RHIC) in BNL, New 
York and the upcoming Large Hadron Collider(LHC) in 
CERN, Geneva. Theoretically, it is expected that the 
spontaneously broken chiral symmetries at low temper- 
atures are restored at high temperatures. While lat- 
tice methods have been very effective in predicting the 
transition temperature as well as the nature of the the 
strongly interacting matter near the transition tempera- 
ture, fermions with exact chiral symmetry on lattice are 
important for such a study of the chiral phase transi- 
tion. Kaplan [l| proposed to define fermions with ex- 
act chiral symmetry on a five-dimensional (5-D) lattice 
with a mass term M in the form a step function(domain 
wall) and with an infinite extent along the fifth dimen- 
sion. The massless 4-D fermions are obtained as local- 
ized on the wall, and are hence known as the domain wall 
fermions. On a finite lattice needed for numerical simula- 
tions, however, fermions of both chiralities exist with an 
exponentially small overlap between the respective chiral 
states 0. Currently, the most popularly used fermions 
in QCD simulations at finite temperatures/densities are 
the staggered fermions which have only a remnant chiral 
symmetry on the lattice. Moreover, they explicitly break 
spin and flavor symmetries. The full chiral symmetry for 
these fermions is recovered only in the continuum limit, 
i.e., in the limit of vanishing lattice spacing. In spite of 
the (exponentially small in N^, the number of sites in 
the fifth dimension) chiral violation on the lattice, the 
domain wall fermions are more promising than the stag- 
gered fermions due to their exact flavor and spin sym- 
metry on the lattice. On the other hand these are more 
expensive to simulate as the computational cost increases 
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linearly with A''5 . One has to optimize A''5 and M for full 
QCD simulations. In order to gain insights on ways to 
minimize the lattice cut-off effects we study various ther- 
modynamic quantities of free domain wall fermions as a 
function of M and with an aim to optimize the irrel- 
evant lattice parameters for faster convergence to their 
continuum values. We find that by adjusting the domain 
wall height M in the range 1.45 — 1.55 rather than the fre- 
quently used choice of M = 1.0, a faster convergence to 
the continuum results for both finite and infinite values 
of A^5 is achieved. However, the cut-off effects are seen 
to be quite large on small lattices with temporal extent 
of 4-6 where most of the current QCD simulations are 
being done. We therefore examine modifications of the 
domain wall, as well as the overlap kernel to minimize 
such corrections for small lattice sizes. 

The plan of the paper is as follows: In section II, we 
compute the energy density of free domain wall quarks 
on the lattice analytically and verify that it yields the 
correct continuum limit. In section III, the same quan- 
tity is computed numerically and the various lattice pa- 
rameters for which the convergence to the continuum is 
fastest are estimated. In section IV, we repeat the calcu- 
lations of energy density in the presence of chemical po- 
tential and susceptibility and confirm that this optimum 
M-range does not shift significantly. In section V, we pro- 
pose a method of reducing the lattice cut-off corrections 
to thermodynamic quantities on small lattice sizes, com- 
puted using both the chiral fermions, namely the domain 
wall at infinite and the overlap fermions. This helps 
in faster convergence to the continuum results even for 
M = 1.0. 



II. ENERGY DENSITY OF DOMAIN WALL 
FERMIONS 

The domain wall fcrmionsil] in the continuum are de- 
fined on a 5D space-time with the mass term in the fifth 
dimension in form of a domain wall (j){M) = Mtanh(s). 
This helps in localizing a fermion of definite chirality on 
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the domain wall. The domain wall operator in the con- 
tinuum is given as, 

4 

£'dw = 5]7m^m + 7555 + </'(M). (1) 

The massless fermion modes in 4-D are obtained when 
the following conditions are simultaneously satisfied. 

4 



^ It^d^ijj = 
(75^5 + 0(M)),A = 0. 



(2) 



It was shown that only one normalizable solution exist, 
bounded to the wall at s = where the 4>{M) changes 
abruptly. The corresponding analog of the domain wall 
term on the lattice is of the form 



<I){M) = Afe(s) 



(3) 



On the lattice we do not get a single massless mode by 
discretizing Eq. ([T]). This is because the lattice regulator 
is anomaly free, so massless fermions of both handed- 
ness exist on the lattice. A Wilson term is needed to 
spatially separate the left and the right handed fermions 
in the 5th dimension by localizing them on the domain 
wall and the anti-domain wall respectively which are sep- 
arated from each other by the lattice extent in the fifth 
dimension N<^. To obtain thermodynamical quantities of 
free fermions with exact chiral symmetry on the lattice 
in 4-D, we need to divide out contribution of the heavy 
fermion modes which exist in the fifth dimension. This 
is done by subtracting a pseudo- fermion action [s^ from 
the standard 5-D action. Following Shamir^], the do- 
main wall fermion action on a x Nt x anisotropic 
lattice with lattice spacings of a, and in the three 
spatial, the temporal and the fifth dimension respectively 
can be written as, 

Sdw = - X! i'jx, s)Ddw{x, s; x', s', fi, m)i>{x , s') 

s,s' — 1 x^x' 



domain wall fermions we do not have a good prescription 
for obtaining the conserved number density. Following 
Bloch and Wetting [J], we incorporate the chemical po- 
tential in Dw but in a general form using the functions 
K and L Q defined below. These multiply the 1 ± 74 
factors in the Wilson-Dirac operator leading to, 



Dw{x,x',fi)= [3+ — -M]S^., 



a I . '^ 1 +74 
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+ L(A)C/4(a:)i^<5,,,,_4 



(6) 



In this paper we consider the non-interacting fermions, 
i.e., U^{x) = 1. Introducing R and 9 by 

^(^l±^^i?cosh. ^i^l^.i?sinh., 
2 2 

(7) 

the free Wilson-Dirac operator in Eq. ^ can be diago- 
nalized in the momentum space in terms of the functions, 

hj = sinapj , = i?sin(a4P4 — iO) , (8) 

3 

/jg = M — — cosapj) (1 — i?cos(a4P4 — z6')). 

a4 



such that 



(9) 



i=l 



To study thermodynamics of fermions one has to neces- 
sarily take anti-periodic boundary conditions along the 
temporal direction. Assuming periodic boundary condi- 
tions along the spatial directions we obtain 



2n.,-7r 



apj 



■^,rij=0,..,(A^-l),j = l,2,3and 



(^^Dw{x,x') + 1^ Ss,s' 

s.s' — l x,x' 

- (P-<5,^.+l + P+Ss.^s-i)5^.x'] s'), (4) 

with the boundary conditions 

F-Vatj+i = -rriqP^^pi, P+Vo = -rfiqP+^pN^ (5) 

where P± = are the chiral projectors and rriq is the 
bare quark mass in lattice units. Dw is the Wilson-Dirac 
operator defined on a 4-D lattice. The volume of the sys- 
tem is V = N^a^ and T = l/{NTa4) is its temperature. 
The chemical potential iia^ = /t is usually introduced as a 
Lagrange multiplier corresponding to the conserved num- 
ber density in the expression for the Lagrangian. For the 



ap4 — LU„ 



N ' ' 

(2n+ l)7r 



,n = 0,..,(A^T-l) (10) 



It is to be noted that M, the height of the domain wall 
on the lattice, is bound to < Af < 2 to simulate one 
flavor quark on the lattice. To suppress the heavy mode 
contributions and recover a single chiral fermion, pseudo- 
fermion fields are introduced which have the same action 



but with 



1 i.e with anti-periodic boundary con- 



dition in the fifth dimension,^]. The fifth dimensional 
degrees of freedom can be inte^ated out to yield an ef- 
fective domain wall operator 6, 7] 

___ ^ 1 + + (1 _ m,)75 Y^^, (11) 
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where the transfer matrix T is 

T=il + ^jr^DwP+)-\l-—l5DwP-)- (12) 
a a 

Since T can be shown to be Hermitian for jj, ~ 0, and 
therefore has real eigenvalues, has only positive 

eigenvalues for even Nr^. Introducing [7] a notation |T|, 



the function 



1+T™5 



in the domain wall operator can be 



expressed in the form of a tanh function as in Eq. (|13p . 



= 1 + mg - (1 - m,)75tanh(-^ In |r|). (13) 



The above derivation of the effective domain wall opera- 
tor assumes that 1 + does not have any zero eigen- 
values. For if it does, then the contribution of the heavy 
modes is zero. If A be an eigenvalue of T, then this as- 
sumption requires that 



InA ^ i 



.{2n + l)n 



(14) 



This is clearly true for /i = for even the interacting 
fermions where T is Hermitian and thus any A is real. 
However, once chemical potential is introduced in the 
Wilson-Dirac operator, as above, Dw and T are not Her- 
mitian any longer for the free fermions themselves, leav- 
ing open the possibility that this condition will not be 
met. 

It is easy to see that three distinct limits are of interest 
in which we should compute the various thermodynamic 
quantities for massless domain wall operator. These are 
as follows: 

1. —^ cxD, as — finite, where one obtains exact chi- 
ral fermions for rfiq — 0, 

2. A^5 — > GO, as — > 0, such that L5 = N^a^, leading to 
an approximate form for the overlap fermions Q, 
and 



different from the overlap case. For this type of Ginsparg 
Wilson fermion too the introduction of chemical poten- 
tial necessarily leads to chiral symmetry breaking [il[ on 
the lattice because the action in presence of /t is not in- 
variant under the chiral transformations Tdj on lattice. 
Like in the case of the overlap fermions, chiral symme- 
try is exactly realized for these wall fermions only in the 
absence of chemical potential. 

The energy density e of the domain wall fermions in 
the chiral limit is evaluated from the temperature par- 
tial derivative of the partition function, Z — T)et{D'j^^,). 
This is equivalent to taking a partial derivative with re- 
spect to a4 on a lattice of fixed size Nt- The energy 
density, 



1 



d 



In Z 



(16) 



can be evaluated analytically in terms of the quantities 
q, s,i and s ,t defined below in Eq. (|17p . where the 
dash denotes the a4-derivative of the respective quanti- 
ties. Defining 



t 



5 ' 



i=l 

s\/s2-4/i5 + 4 , q = s'^-2h5 + 2, 
dh4 a 

— — = —Hsm(a4P4 — lO) , 
aa4 aj 



daA 



r(l — i?cos(a4P4 — id)) 



s 

s t s^(ss — 27) 
s t 



(17) 



one has 



3. — finite, = finite, corresponding to the form 
of domain wall operator directly relevant for prac- 
tical simulations on the lattice. 



oo,as = finite 



In this limit, the tanh function in Eq. (|13p becomes 
sign function and the resultant effective domain wall op- 
erator is given as 



1 



2t 4/i57 + 4ss (1 + 5^) 



2hl + 2s2 + s4 



N^Nt ^ \ t 

Pj.n \ 

+2ss't + s^t - 475^ - 8/15SS' - 27t - 2h^t 



-sH - 4/i5s2 _ 2h^t 



1 



J2FiR,^n,p). 



Pj,n 



(18) 



D 



eff _ 
DW ~ 



l + rhq- {1- mg)75e(ln |r|) 



(15) 



For rhq — 0, this form of the domain wall operator sat- 
isfies the Ginsparg- Wilson relation Q, as shown in the 
appendix [XI Indeed, it is just like the overlap operator, 
but with a different argument of the sign function. The fi- 
nite size corrections to various thermodynamic quantities 
computed with this lattice operator are expected to be 



Setting a = a4 after evaluating the a4-derivatives, the 
summation over the discrete Matsubara frequencies can 
be evaluated analytically by the standard contour inte- 
gral technique or numerically by explicitly summing over 
them and the momenta pj . For the former, we need to de- 
termine the singularities of the summand F in Eq. (jlSp . 
We outline below briefiy the results one obtains for the 
zero and finite temperature cases. 



4 




FIG. 1: Contour chosen for evaluating the energy density for 
nonzero value of chemical potential at zero temperature. The 
thick line indicates the Matsubara frequencies while the filled 
circles denote the poles of F{R,uj). 



T = 0, fi ^ : In order to obtain a general condi- 
tion for eliminating the spurious /t^-divergences, we first 
calculate the energy density at zero temperature in the 
limit Nt — > oo at finite a. The frequency sum 1 /Nt 
in Eq. p8|) gets replaced by the integral ^ J^^ dw in 
this limit. Subtracting the vacuum contribution corre- 
sponding to /t = 0, i.e. R = 1,9 = 0, the energy density 
at zero temperature is given by 



Pi 



F{R,uj-i6)duj- / F{uj)du 

-TT 'J —TT 



For brevity, we suppress from now on the explicit pj- 
dependence of the function F although we retain the 
overall sign to remind of it. Choosing the contour shown 
in Figure [U the expression above can be evaluated in the 
complex oj-plane as 



ea 



-y 



27ri^Res F(R,uj^) 

T 

F{R,uj)duj 



F{R,uj)dw 



TT — tU 



F{R,uj)dLo 



F{uj)du 



(20) 



The second and fourth terms cancel since F satisfies 
F(i?, TT + irj) = F{R, — TT + ii]). Hence, we obtain 



1 



ea 



E 



Pj 

F{R,uj)duj 



F{uj)duj 



(21) 



where -ii?i is the residue of the function F{R, uj) at the 
pole —i smh~^ {y/J / R) . It is clear from Eq. ([2T|) that 
the condition R = 1 cancels the integrals, yielding the 
canonical Fermi surface form of the energy density. For 
R 1, there will in general be violations of the Fermi 
surface on the lattice. Moreover, in the continuum limit 
a — *■ 0, one will in general have the /^^-divergences for 
i? 7^ 1 in the energy density. The condition to obtain the 
correct continuum values of e = /i^/47r^ turns out to be 
K{fi) — L{jj,) = 2jl + 0{ij.^). That this effective domain 
wall fermion satisfy the same condition as the overlap [llj 
suggest that such condition may be generically true for 
Ginsparg- Wilson fermions. Also that one obtains iden- 
tical condition in the staggered caseQ suggests that the 
behavior near the continuum limit dictates this condi- 
tion. Note also that the form used by Bloch and Wettig 
namely, exp(±/i) for K, L, also satisfies the condition 
R^K -1=1. 

T ^ : In order to choose the appropriate contour 
in the T ^ case, note that the function F{R,uj,p) 
has poles at cos^^{y^d'^ — g) = ±ismh^^ y/f. These 
turn out to contain the physical poles in the contin- 
uum limit. However, there are additional(unphysical) 
poles at cos~^(— -^/d^ — g) — ±7r ± i sinh""'^ ^//i ±7r ± 
icosh"^ ^(zbicosh^^ ^) for ^ > 0(< 0) and at 

±icosh~^a;i where = It^^n-f ■ The definitions of the 

2(3-2) 

quantities d,f,g are 



9 
b 

f 

d 



M-A + b, with 
cos(api) + cos{ap2) 
hl + hl + hi 



cos(ap3) 



= 4 + (Af - 4)2 + 2(M - 4)6 - 



with 



i<i<4 



2 cos{api) cos(apj) 



(22) 



Note the obvious similarity with the case of the overlap 
fermions IjJ , where the same set of poles contribute as 
well. Unlike in the overlap case, however, there are no 
branch cuts present in this case. The contour chosen 
for evaluating the frequency sum shown in Figure [H is 
thus different from that chosen for overlap fermions. The 
residue of the pole enclosed by the contour for F comes 
out to be, 



^ Vf VT+7-1 
VTo+T) 



G(M). 



(23) 



with the first term yielding the continuum value of the 
energy density in the limit of vanishing lattice spacing 
a. The energy density expression comes after performing 
the contour integral comes out to be. 



ea 



-y 



Pi 



1 



^Nt sinh ^77-1-1 



V7 I X/TT7-1 
\/i+7 v//(i + /) 

+ £3 + £4 , 



G(M) 



(24) 



5 




FIG. 2: Contour chosen for evaluating the energy density 
at finite temperature. The crosses indicate the Matsubara 
frequencies while the filled circles denote the poles of F{p, uj). 



operator, the energy density can be evaluated[l3| as, 



Pj ,n 



4sinh[^]((-/i4/i5a + h^-f) cosh[^] 



sN^Nj 



+ {hih^a + [hi + s^)-i + 2h^sH) cosh(^) - 2ssinh[4^ 

[h^ + (s2 + /i2)cosh[2sL5] 

{hit + /t57 + [h'^t + hja + 2h5{hzt + 7)) cosh[sL5])) 
—2h^ssmh\^sL^]) 



(26) 



where a and 7 are the same as defined previously and t 
is now defined as, 



t = 



(- sin^ api + /i57)(- tanh + -%^sech^^) 



tanh 



(27) 

It was checked that the overlap energy density is obtained 
back when L5 — > cx3. We use the expression above for our 
numerical work presented in section III. 



which again turns out to be similar to the overlap case. 
Due to a different functional form of F and a different 
choice of contour, the corresponding lattice correction 
terms £3, £4 which are the line integrals of F along lines 
3,4 in the Figure^ are different, leading to different finite 
size corrections. In the continuum limit the unphysical 
poles are pushed to infinity and the values of £3, £4 vanish, 
leaving only the contribution of the physical poles to the 
energy density: In the square bracket, only first term 
gives the usual continuum expression with the other term 
vanishing as a — > 0. The same treatment goes through in 
presence of /t only the contour has to be shifted along the 
imaginary lo plane by an amount dependent on (1 with the 
position of the poles in the complex w-plane remaining 
unchanged. 



C. Finite A^s and 05 

While performing Monte Carlo simulations with do- 
main wall fermions one needs to work on lattices with 
finite number of sites in the fifth dimension. For finite 
iVs, the chiral symmetry is broken and it is important 
to ascertain the dependence of the correction terms with 
iVg. Evaluating the matrix tanh(iV5/2 In |T|) in Eq. (fT^ 
various thermodynamic quantities of free domain wall 
fermions on the lattice can be evaluated. The energy 
density in the massless limit then is 



ea 



Pj, 

tu + ut 



22^5 +1 +2^=u 
xq — (<Z — 2)x \ 



B. N. 



finite 



tu — (q — 2)x 



(28) 



In the case when the lattice spacing in the fifth direc- 
tion a5 — > and the number of sites — > 00 such that 
L5 — N^ac, is finite, the effective domain wall operator 
reduces to 



Ddw = (1 + nig) + (1 - mg)jH(inh{—Y''Dw) 



(25) 



where the quantities u and x are functions of h's defined 
in Eq. ^ defined as, 



t 



h5-l 
t-q 



1 



t + q 

1-/15 

t + q 
1-h. 



N5 



(29) 



Starting from the above expression we recover the overlap 
operator when L5 ^ 00. With this effective domain wall 



The partial derivatives of the above variables are repre- 
sented as the same variables with a dash, and are func- 



6 



tions of h's, a and 7. 
q — 2ss 
u = N5 
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t-q 


) 


<' - 


- q 


7(t 


-9) 




1 


115 - 


— i 


("5 - 


1 ^2 


t + q ^ 


-/V5-I 

) 


<' 4 


-f + 


7(i 


+ 9) 




/ 


1 - 




(1- 




t^q ^ 




t - 


- q 


7(t 


-q) 


h5-ly 




h - 


- 1 




-1)2 


t + q'' 




t 4 






+ 9) 


1-/15, 




1 - 


7 + 

/l5 


(1- 





Again, we shall use these expressions for obtaining the 
numerical results presented below where we also show the 
results for quark number susceptibility. The same set of 
formulae remain valid for calculation of the susceptibility 
except for the fact that and 7^ then are the derivatives 
with respect to p, and are defined as, 

dh^ ia 
ctf^ = -7r:r = — cos(a4P4 - »m) , 
ojl 04 

7p = — ^ = —ihn{ioY number density) 




djl 



III. NUMERICAL RESULTS FOR A = 



FIG. 3: The C, dependence of the energy density of domain 
wall fermions for M = 1.50, as = 1 and in the limit A'5 — > 00. 



TABLE L e/esB values for different M for C = 4 



Nt 


M=1.0 


1.40 


1.45 


1.50 


1.55 


4 


0.909 


1.240 


1.286 


1.333 


1.384 


6 


1.308 


1.413 


1.426 


1.444 


1.469 


8 


1.317 


1.221 


1.197 


1.175 


1.159 


10 


1.236 


1.090 


1.051 


1.009 


0.966 


12 


1.168 


1.048 


1.011 


0.966 


0.915 


14 


1.123 


1.043 


1.014 


0.977 


0.930 


16 


1.092 


1.046 


1.026 


0.998 


0.960 



A. A^5 = 00, as = 1 

The goal of our numerical study is to find the opti- 
mum range of M for which the finite lattice spacing cor- 
rections are minimum and compare it with that for the 
Dirac-Neuberger case [lH. We do this in the chiral limit 
and set rhq = 0. The lattice energy density given by 
Eq. psp was computed numerically by summing over the 
momenta along the spatial and temporal directions. The 
zero temperature part of the energy density was deter- 
mined in the limit Nt — > 00 on a lattice with a very large 
spatial extent N by numerically evaluating the ap4 = lj 
integral. Holding the physical volume constant in units 
of T by keeping V^^^T = N/Nt = C fixed, we define 
the continuum limit by Nt — > 00. The thermodynamic 
limit is then achieved in the limit of large C- We first 
determine the acceptable range of C by looking for (- 
independence. The e obtained by subtracting the zero 
temperature part from the lattice energy density expres- 
sion was normalized by its continuum value ess- Figure 
[3] displays the ratio e/esB as a function of Nt for differ- 
ent values of C at a fixed M — 1.50. One notices that for 
C > 3 the energy density plots lie on top of each other, 
suggesting the thermodynamic limit to have reached by 
C = 4 — 5. In order to highlight the deviations in the 
continuum limit, the same ratio is exhibited for different 
M values for C = 4 in Figure [Has a function of l/N"^ and 



listed in Table|T]for a range of Nt likely to be used in sim- 
ulations. We choose to define the optimum range of M as 
the values of M for which the thermodynamic quantities 
are within 3% of the continuum values for the smallest 
possible Nt- One sees from both the Figure [4] and the 
Table H] that the order 1 /N^ corrections are minimum 
for M between 1.45-1.50 and Nt > 12. The correction 
terms for Af = 1 are linear in 1/N^ for Nt > 10 and 
are about 20% of the continuum value even for Nt — 12. 
This is similar to that reported earlier for the overlap 
fermions [ll|. Though the continuum extrapolation with 
M = 1 is easier due to the linear functional form, it is 
computationally expensive, needing simulations at more 
values of Nt, each greater than 10. Thus M =1.45-1.50 
seems to be an optimum range for lattice simulation of 
the energy density of domain wall fermions. We have 
also varied the lattice spacing along the fifth dimension 
05 to find out how the cut-off dependent terms change 
with it. The correction terms to the energy density for 
a5 = 0.5 at small lattice sizes Nt < 10 are indeed larger 
than that for a<^ — I for the above mentioned optimum 
range but for iV^ > 12 such terms are again within 2- 
3% of the Stefan Boltzmann value. The optimum M 
range for which the lattice artifacts are minimum shifts 
to 1.50-1.60. Thus there is a marginal dependence on 05 
for Nt > 10. Reducing 05 further does not increase the 
range much as we demonstrate in the plot for 05 ^ in 
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FIG. 4: The variation of energy density of domain wall 
fermions with M in the limit Afs — > oo and as = 1. 
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FIG. 6: The variation of energy density on lattice with 1/A^|. 
for domain wall fermions for L5 = 14 and different M. 
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FIG. 5: The variation of energy density on lattice with Nt 
for domain wall fermions for different Ls, as shown by the 
respective labels, and M — 1.55. 



Figure [6l 



B. Ai's — > 00 , as — > , L5 = A'"5a5 — finite 

Next we investigated the limit ^ oo , 05 — > such 
that L5 = A^5a5 = finite in order to estimate numerically 
the value of L5 for which we recover the overlap energy 
density starting from Eq. ((26)) . As can be observed from 
Figureini Ls-independent results are obtained for L5 > 14 
for M — 1.55. This was also the case for a range of AI 
around this value. For L5 < 10 the convergence towards 
the esB value was seen to be very slow for all M and 
we find that the continuum value is not reached even for 
lattice size as large as Nt = 32. Figure [S] displays the 
results as a function of l/N^ for L5 — 14 and various 
values of M indicated on it. The deviations from the 
continuum for such L5 are less than 3 % for the range 
of M between 1.50-1.60, in agreement with the overlap 
results [ni- 



TABLE II: e/ess values for different M for C = 4, N5 = 18 



Nt 


M=1.0 


1.40 


1.45 


1.50 


1.55 


4 


0.909 


1.240 


1.285 


1.333 


1.383 


6 


1.308 


1.413 


1.425 


1.443 


1.467 


8 


1.317 


1.221 


1.197 


1.174 


1.156 


10 


1.237 


1.090 


1.052 


1.009 


0.965 


12 


1.169 


1.049 


1.013 


0.968 


0.917 


14 


1.123 


1.045 


1.017 


0.980 


0.934 


16 


1.093 


1.048 


1.029 


1.002 


0.966 



C. Finite and 05 = 1 

The case of finite with as = 1 is clearly of 
most interest for practical simulations with dynamical 
fermions. Earlier numerical studies for free domain wall 
fermions [Tl, [l3| employed M = 1.0 and found some- 
what slow convergence of various thermodynamic quanti- 
ties towards their continuum values. We intend to check 
whether tuning the value of M results in a faster conver- 
gence. For that purpose we have computed the energy 
density expression for finite N5 and 05 = 1 in Eq. ([28)l by 
summing over all the discrete momenta. We display those 
results for e/esB in Figure [T] The upper panel shows the 
results for a series of and a fixed M — 1.5. The results 
are seen to become TVs-independent by — 18, making 
it an optimum choice for obtaining continuum results on 
the lattice. The lower panel shows the M- variation for 
= 18. Table |TT] provides the values for lattices with 
reasonable iVr-extent. The general trend is clearly the 
same as above with M =1.45-1.50 emerging as the range 
for which the Stefan-Boltzmann limit is reached to within 
3-4% for Nt > 10 (Table |TT|. Interestingly, N5 = 18 
seems to mimic the —^ 00 limit quantitatively rather 
well as can be seen by comparing the Tables |T] and [TTl 
Consequently, the same optimum range of M is obtained 
for both. 
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FIG. 7: The variation of energy density on lattice with 

for domain waU fermions at different for M = 1.50 (upper 

panel) and at = 18 for different (lower panel). 



FIG. 8: The energy density of domain wall fermions in pres- 
ence of /i for different M for fi = 0.5/Nt and infinite A^s (upper 
panel) and for A^s — 18 (lower panel). 



IV. NUMERICAL RESULTS FOR fi ^ 

It should be noted that in this case T is no longer 
Hermitian but as long as the condition given in the Eq. 
p4)) is satisfied the effective operator in Eq. (flSjl is well 
defined. We shall restrict the range of fi to ensure that 
it is so. We choose K and L to be e^'^ respectively in 
our numerical computations as suggested in [4| . Our aim 
again is to find the optimum M for which the continuum 
results are obtained with least computational effort, and 
compare it with our the range obtained from the energy 
density above. We consider two observables here. One 
is the change in the energy density due to nonzero fi : 
A.e{fi, T) = e(/i, T) — e(0, T). In the continuum limit this 
is 



Ae(^,T) 



j£_ 

T4 Att^T^ ^ 2T2 



(30) 



Another observable we studied was the quark number 
susceptibility at /I = 0. It is defined for any jl by, 



1 



X = 



IndetD 
d{j? 



(31) 



and in the continuum is given by, 

,2 j,2 



(32) 



We will focus on x(0) here due to its importance in the 
applications to the heavy ion collisions. 

We estimated numerically Ae(/i,T) for fi/T — jlNx 
fixed at 0.5. The upper and lower panels of the Figure 
[5] display our results for this observable in the units of 
T'^ for iVs = oo and 18 respectively for the M values 
indicated. The horizontal line in each case shows the ex- 
pected result in the continuum limit from Eq. pop . From 
the Figure [5] it is evident that there are no /j,^ diver- 
gences on the lattice, as expected. The deviations from 
the continuum limit are due to the M dependent finite 
size effects. These correction terms are again seen to be 
small for the same optimum range of 1.45 < Af < 1.50 for 
both the cases, as obtained in the zero chemical potential 
case in section III. 

The A^s dependence of the quark number susceptibil- 
ity at /i = is plotted in Figure [TOl It too exhibits a 
convergence to the infinite results for > 16, indi- 
cating that TVs = 18 can again be used safely to approxi- 
mate the infinite N5. Figure [TO] shows the M-dependence 
of the quark number susceptibility at /i = 0. Both the 
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tions, which should however be small for large enough 
temperature where one expects those computations to 
approach the free quark gas results. In all our plots we 
find that for the optimum M range, the deviations from 
the ideal gas results at smaller Nt =4-8 are quite signifi- 
cant but with a relative mild M-dependence for M > 1.4. 
Thus a slightly larger value of M than the optimum range 
we found may not change the finite size effects drastically 
for small Nt- What one does need to be careful about 
though is the extrapolation to the continuum limit. For 
the optimal range of M and Nt > 10, the smallness of 
corrections compared to other errors in the computations 
may make it a less important issue. 



FIG. 9: The quark number susceptibility as a function of 
l/N'}' for values as indicated for AI = 1.5 and C ~ ^■ 



IMPROVEMENT OF THE CHIRAL 
FERMION KERNELS 




FIG. 10: The quark number susceptibility as a function of 
1/Nt for M values as indicated for C = 4 and JV5 = 00 (upper 
panel) and 18 (lower panel). 



Nq — 00 (upper panel) and 18 (lower panel) show small 
deviations from the Stefan-Boltzmann value of 1/3 for 
1.45 < M < 1.55 range and for Nt > 10. Recent 
computations of this susceptibility[lg] for the interacting 
domain wall fermions were performed with Af = 1.8. Of 
course, one expects some shift in M due to the interac- 



In the previous sections we observed that the fermions 
with exact chiral symmetry on the lattice have large 
1/N!^ corrections for small Nt- While we found that 
the continuum limit for various thermodynamic quan- 
tities can be approached faster by choosing the irrele- 
vant parameter M in the range 1.45-1.55, the correction 
terms for Nt =4-6 are about twice that of the Stefan- 
Boltzmann result for such a choice of M too. Here we 
describe our attempts to improve the convergence to the 
continuum results for small Nt and even for A I = 1.0. 
Having the option of the choice of M = 1.0 may be useful 
since it has been noted previously [3, [l3| that the residual 
mass for such a choice of M is zero for a range of N^ at 
the tree level. 



A. Domain wall kernel 

The domain wall operator given in Eq. (jl3p is a 
matrix-function of the Wilson-Dirac operator as in Eq. 
(|6]) . It is clear that its improvement may lead to a better 
domain wall operator, or indeed even a better overlap 
operator, one is looking for. Inspired by the attempts 
to improve the staggered fermions in the so-called Naik- 
action [l^l , we add three- link terms to the Dw as below. 



Dw{x, x',ii) = [3 + — - M ] 5^^^, 



1 - ci7j 



X^X —J 



(33) 
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FIG. 11: The susceptibility of improved and the conventional 
domain wall fermions at M = 1.0 as a function of l/N^. 



Comparing with the Eq. ([S]), it is clear that the modi- 
fication amounts to replacing 7^ by (ci + 03/8)7^. The 
Wilson mass term, added to remove the doublers, is kept 
unchanged. Note that the modified _DvF-operator is still 
75-hermitian for arbitrary real values of the coefficients 
ci and C3. The new domain wall operator can therefore 
be derived in the same way as Eq. (|13p was obtained. 
We fix the coefficients by demanding the dispersion re- 
lation for free fermions on the lattice to be the same as 
in the the continuum up to 0(a*p*). We find that all 
the terms at O(a'^p^) are eliminated for the coefficients 
ci = 9/8, C3 = —1/8. We employ them below for the 
calculation of the thermodynamic quantities. 

The ratio of quark number susceptibilities, x/xsb, 
computed using the above modified domain wall oper- 
ator [l§|, is plotted as a function of as in Figure [TT] 
along with that for the unimproved domain wall operator 
of Eq. 1131). We used M = 1, C = 4, 7V5 = 18 and 05 = 1 
for this computation. One clearly notices that the large 
correction terms 45%) at Nt =6-8 for the usual do- 
main wall operator go down to about 7-8 %. Indeed, the 
size of corrections go down further as Nt increases. Sim- 
ilarly, the the energy density of such improved fermions 
also exhibited smaller, about 15-5%, deviations from the 
continuum for Nt =6-10, as compared to about 30 % in 
the lower panel of Figure [T] 



B. Overlap kernel 

From section II, we know that the overlap operator 
can be derived as a special limiting case of the domain 
wall operator. It would be thus interesting to check how 
the improvement in the Wilson Dirac operator in Eq. 
([55)1 fairs in the overlap case. For that purpose we com- 
pute the quark number susceptibility for non-interacting 
fermions on a N^ x Nt lattice numerically with the cor- 
responding improved overlap operator. The x/xsb does 
have lower 1/N^ corrections for Nt = 6,8 than for the 




0.005 0.01 0.015 0.02 0.025 0.03 

1/n/ 



FIG. 12: The x/xsb as a function of l/N^ for the improved 
and the canonical overlap operators at M = 1.0. 



conventional overlap operator with M = 1 as shown in 
Figure [121 We also observe a faster approach to the con- 
tinuum result with such improved overlap operator than 
with the Neuberger overlap operator even with optimum 
M = 1.55 reported in [ll|. Another advantage is that the 
thermodynamic quantities calculated from this improved 
operator are free from oscillations at odd-even values of 
Nt exhibited JlJ by the usual overlap operator. The im- 
provement in the energy density is marginal up to Nt = 8 
but substantial for Nt = 10 onwards. 



VI. CONCLUSIONS & DISCUSSIONS 

Since the chiral violations vanish exponentially with 
the number of sites N^, in the fifth dimension, the do- 
main wall fermions offer a more practical alternative to 
the overlap fermions and yet have exact flavor and spin 
symmetry. We have computed the energy density and 
susceptibility at zero chemical potential of such fermions 
numerically for both finite and infinite A^s. The chiral 
symmetry is exact in the latter case and a choice M be- 
tween 1.45-1.50 allows faster convergence to the contin- 
uum results. We have also verified analytically that the 
energy density has the correct continuum value in the 
chiral limit. Varying the number of lattice sites in the 
fifth dimension, we have shown that A'^s = 18 is sufficient 
to restore chiral symmetry. 

We found that introducing chemical potential jl in do- 
main wall operator leads to chiral symmetry breaking 
even for infinite A^s. But if we do allow that, there 
exist a large class of functions K(fi) and L(fl), with 
K{jl) ■ L((i) — 1, for which there are no /i-dependent 
divergent terms in the physical observables. From the 
numerical evaluation of the energy density in presence of 
/t , we conclude that the optimum range of M remains 
the same. The lattice cut-off effects are however very 
large for small A^^ = 4-8. By systematically removing 
the dominant correction terms to the continuum value 
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of the chiral fermion operators we have achieved a faster Hence {75, D]jw} = DowlbDow 
convergence to the continuum as well as small 1 /N^ cor- 
rections for small lattice sizes even for M = 1.0. This set 
of optimum parameters is anticipated to produce simi- 
lar results in full QCD simulations with chiral fermions 
though an explicit check needs to be done. 
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APPENDIX A: PROOF FOR GW RELATION 

We show here that the effective domain wall operator 
Eq. (fT5|l satisfies the Ginsparg Wilson relation : 



75DDW + Dowlb = 75(1 - 75e(ln |T|)) 

+ (l-75e(ln|T|))75 

= 2 75 -e(ln|T|)-75e(ln|r|)75 

DdwI^Ddw = (l-75e(ln|r|))75(l-75e(ln|T|)) 

= 75 - 75e(ln |r|)75 - e(ln |T|) + 75e2(ln |r|) 



Since, 
e(ln|r|) 



InlTl 



^{\n\T\)\n\T\ 



e^lnT) = — ' ' ' ' - = 1 Al 
(^ln|T|ln|r|)2 
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